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Abstract 

We investigate the aging behavior of lattice-gas models with constrained 

dynamics in which particle exchange with a reservoir is allowed. Such models 

provide a particularly simple interpretation of aging phenomena as a slow 

approach to criticality. They appear as the simplest three dimensional models 

exhibiting a glassy behavior similar to that of mean field (low temperature 

mode-coupling) models. 
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As the temperature of a supercooled liquid is reduced, its dynamics becomes slower and 
slower and its viscosity increases. The glassy transition is usually located at the temperature 
Tg where the viscosity first exceeds the conventional value of 10^"^ Poise At lower temper- 
atures, the glass is usually assumed to be "solid" in the sense that its molecules rattle within 
a "cage" formed by their neighbors and do not move away from it, at least within experi- 
mental time scales. However the glassy state is not an equilibrium one since the evolution of 
the system does not stop altoghether, but rather keeps evolving at a slower and slower pace 
as the time tw elapsed since its quench increases. This is the origin of the striking aging 
effects observed in glasses |0J^. Aging corresponds to the property that, while one-time 
quantities like the average energy, volume, etc., appear to be invariant under time transla- 
tions, two-time correlations and responses exhibit a non trivial dependence on both of their 
arguments 0. Similar properties have also been observed in spin-glasses 0. These aging 
properties have been thoroughly investigated within spin-glass models |^ which are closely 
related to the mode-coupling theories of structural glasses [Q. Aging appears provided one 
takes the infinite size limit before taking the infinite time limit In this situation, the 



Boltzmann-Gibbs equilibrium distribution is never reached. Conventional mode-coupling 
theories, which are based on the hypothesis that the initial distribution is the equilibrium 
one, are therefore inconsistent in this limit below Tg. 

The spin-glass models proposed as a description of structural glasses lack, however, a 
transparent physical interpretation in terms of particles and involve a complex, random 
hamiltonian which is hard to justify as a description of a fiuid. Their justification is rather 
a posteriori in the sense that the phenomena they exhibit recall the behavior of structural 
glasses. On the other hand a class of very simple kinetic models have been introduced to 
describe the slowing down of the dynamics [|ri|,|T^. These models are defined by kinetic 
rules involving a selection of the possible configuration changes ( "moves" ) and are therefore 
called models with constrained dynamics. The kinetic rules satisfy detailed balance and 
are compatible with a Boltzmann-Gibbs equilibrium distribution involving a hamiltonian, 
usually chosen to be a trivial one. The constraints are alone responsible for the slowing down 



of the dynamics because, near any allowed configuration, there are only few configurations 
which satisfy them. 

Some out of equilibrium properties in models of this kind have been already discussed 
in [|13|. We shall in this paper refer to the kinetic lattice-gas model studied by Kob and 
Andersen (KA) |T^. The system consists of particles in a cubic lattice of side L, {V = L^) 
with periodic boundary conditions. There can be at most one particle per site. At each 
time step a particle and one of its neighbouring sites are chosen at random. The particle 
is moved if the three following conditions are all met: (i) the neighbouring site is empty; 
(ii) the particle has less than four nearest neighbours; (iii) the particle will have less than 
four nearest neighbours after it is moved. The rule is symmetric in time, detailed balance 
is satisfied and the allowed configurations have equal probability in equilibrium. In this 
model the dynamics becomes slower and slower as the particle density p increases. Above a 
critical value of the density, pc, particle diffusion stops completely. These investigations are 
carried on initializing the system in an "equilibrium" configuration: therefore the approach 
cannot be extended below the glass transition temperature without becoming vulnerable to 
the same objection voiced above against the mode coupling theories: how can one reach this 
"equilibrium" configuration? 

We show in this paper, by numerical simulation, that it is possible to study aging in 
these models, if they are provided with a mechanism allowing one to realize a process analog 
to a quench. We consider the three-dimensional KA model, but we allow the lattice gas to 
exchange particles on a single two-dimensional layer with a reservoir characterized by the 
chemical potential /i [1^ . For each value of /x it is trivial to evaluate the equilibrium value of 
the density Peq{p)- If Peq{p) < Pc, the system rapidly reaches the equilibrium state. There is 
therefore a critical value, pc, of p, defined by PeqiPc) = Pc- A quenching process corresponds 
to performing a jump of p from below to above pc- Therefore, p plays a role analogous to 
the inverse temperature in mode-coupling theories. We then observe that p never exceeds 
Pc, but rather approaches it like a power law in time, and that the square displacement and 
the self-correlation function of the particle violate time-translation invariance, since they 
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explicitly depend on the waiting time, t^,, after the quench. 

One can also perform "annealing" (compression) experiments, where fi is increased at a 
fixed rate. In close resemblance with the behavior of real glasses, we observe that the limit 
density depends on the compression rate, and tends to pc as the rate decreases to zero. We 
can then evaluate, by integration, the entropy variation of the reservoir. The curves remain 
consistently above the equilibrium curve S = Seq{^), and seem to reach a critical value, Sc, 
as the compression rate decreases to zero. This is reminiscent of the Kauzmann paradox, 
which is usually considered as an argument for the existence of a thermodynamical phase 
transition related to the glass one. However in this case, we know that there is no such 
transition. 

Hysteresis effects also appear when the reservoir's chemical potential is varied cyclically. 

As mentioned in the Introduction, KA studied the diffusion and the decay of the corre- 
lation functions evolving at fixed number of particles starting from a random (equilibrium) 
configuration. They found that correlations decayed to zero with a characteristic time r(p) 
that appears to diverge at ~ 0.881, where also the diffusivity constant goes to zero. For 
densities above pc, the correlations did not decay at all. 

We modify the KA model by creating and destroying particles on a single slice, say 
{x,y,z) = {l,y,z) (the "surface") with the following Montecarlo rule: We choose a site on 
that slice at random. If it is empty we add a particle, if it is occupied we remove the particle 
with probability e~^. We alternate such sweeps of creation/destruction with the ordinary 
diffusion sweeps. Since the creation/destruction mechanism is very fast, the surface itself is 
always in equilibrium. 

It is easy to calculate the thermodynamical equilibrium quantities (we set the temper- 
ature ksT = 1 throughout): 1/p = 1 + e^^, where p is the density. The entropy per unit 
volume is given by 5* = —pin p — (1 — p) ln(l — p). 

We first consider compression experiments in which 1/p (a quantity analogous to tem- 
perature) is slowly lowered from a high value to zero. In Fig. |T]a we report the results of the 
specific volume v = 1/p versus 1/p for several annealing speeds. The smooth curve is the 
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equilibrium state equation and the horizontal dashed line is the critical specific volume of 
KA Vc = ^-1 pc = 1/0.881. In Fig. ^d, we show the change in entropy of the reservoir, given 
by 5'(yu) = S'(yUj) — fidp. The smooth curve is the equilibrium curve S vs. 1/ p, and the 
horizontal dashed line is the entropy at the critical density. These curves look quite similar 
to the energy vs. temperature and entropy versus temperature curves in the annealing of 
glasses. As predicted, the system falls off equilibrium around pc- We have checked that 
the bulk of the sample does not exhibit density inhomogeneities except very close to the 
"surface" . 

We now turn to the behaviour of the system after a quench to the subcritical value 
1/fi < l/pc- Here we concentrate on the behaviour of a large system at long but finite 
times. In Fig. |^ we plot the specific volume minus its asymptotic value versus time after 
a quench to the subcritical value l/p = 1/2.2. The data can be fitted by a power law of 
the form v{t) — Voo ~ where a = 0.57 ± 0.02. The asymptotic value of the density 
Poo = l/'i^oo = 0.861 is slightly lower than the KA value pc = 0.881, suggesting that the 
critical density might have a weak dependence on the way the particles entered the bulk 
(although it might also be a finite-size effect). 

The mean square displacement R'^(t + tw,tw) of a particle between time and time 
t + ty^ must be defined with some care, since the particles may leave or enter the system. We 
define it by averaging only over the particles which are present at both times. We find that 
diffusion indeed slows down as the waiting time increases. However, the time t{R^) needed 
to reach a given square distance is roughly proportional to only for long times. In fact. 



we find here a phenomenon already encountered by Kob and Barrat in Lennard- Jones 
systems: one has to consider an effective waiting time that takes into account the relaxation 
of the system that happened already before the quench. Hence, one should define an effective 
waiting time f^-^ = + tq, where tq is the relaxation time characteristic of the equilibrium 
situation before the quench. Fig. ^ shows that, if one plots the squared displacement vs. 
t/tlJ-l' , the curves lie roughly on top of each other. The deviations are probably due to finite 
size effects, e.g. to particles that escape from the sample. The plot exhibits the waiting-time 



dependence of the diffusion constant and is consistent with a simple argument suggesting 
that diffusion is logarithmic in the presence of aging 



One can also compute easily the effective potential of Franz and Parisi |T6|. This potential 
is a purely static quantity, and therefore is independent of the dynamical constraints. Indeed, 
we have not found any signature for the dynamical transition in the effective potential. 
However the "number of neighbours" constraint is an idealization of microscopical (repulsive) 
energy terms, and as such it should be seen by the effective potential. It would be interesting 
to clarify this point further. 

Since the size of the system considered here is finite, equilibrium must eventually be 
reached, provided that any two allowed configurations can be connected by a path of allowed 
moves. It is possible to convince oneself that the kinetic rules allow an initially empty lattice 
to be progressively filled in, leaving only 0(1/ L) empty sites per unit volume, and that it 
is possible to find a path connecting almost any two allowed configurations, if necessary by 
letting the particles escape one by one by the way they got in. Therefore, for any value of 
the equilibrium density will eventually be reached, but with times which diverge fast 
as L ^ oo. We have verified that a system of size 5^ equilibrates after a time of order 10^ 
MC sweeps. What this argument implies is that after a quench there is first the power-law 
approach to a critical density pc, followed by a much slower (size-dependent) increase in 
density to the equilibrium density. 

We have thus shown that the out of equilibrium dynamics of kinetic lattice-gas models 
reproduce the phenomenology of more sophisticated glass models if one endows them with 
a surface mechanism of particle exchange with a reservoir. Equivalently, we could have 
coupled the system to a piston and used the pressure as an external parameter instead of 
the reservoir's chemical potential. The qualitative behaviour of this model is similar to that 
of mean-field spin-glass models, whose dynamical equations are exactly (in the infinite-size 
limit) the "idealised" mode-coupling equations and their correct low-temperature extension. 
The roles played by the inverse chemical potential and the inverse density (specific vol- 
ume) in the present model are analogous to the ones played by temperature and energy 



respectively, in the mean-field case. When quenched below Tc, mean-field spin-glass models 
never equilibrate. Their energy relaxes to a "threshold" value Ethres{T) which is higher 
than the equilibrium one ||T0|, and can be characterized as the value below which the phase 
space breaks into disconnected components. In our case, the moves allowed by the kinetic 
constraints keep the phase space of the system connected only for p < p^. Therefore the 
density approaches this threshold value, and the corresponding manifold of allowed states 
approaches to the incipient "percolation" cluster characteristic of pc [^. Only for finite 
sizes the system is able to reach values of p larger than pc, just as energy values below the 



threshold can only be reached by finite systems in mean- field spin-glasses [|T0|,|T9|. 

In short: in both cases the system ages because it slowly evolves towards a situation 
in which the available phase-space becomes completely disconnected pO|. The closer the 
system to the "phase-space percolation" threshold level, the slower its dynamics becomes. 
In this sense the system ages because it approaches criticality |^ . 

To summarize, lattice gas models with constrained dynamics, in spite of their simplicity, 
capture the essential features of the aging behavior of three dimensional glasses. 

We are indebted to J.-L. Barrat for suggesting to look at the effective potential and 
to J.-P. Bouchaud, L. Cugliandolo, P. Holdsworth and W. Kob for useful discussions and 
suggestions. M.S. acknowledges the hospitality of the Laboratoire de Physique Theorique of 
Ecole Normale Superieure de Lyon, where part of this work was performed. 
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FIGURES 
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FIG. 1. a: Compression experiment. The smooth curve is the equilibrium state equation; while 

the dashed line is the KA critical value of specific volume, b: 'Entropy' as obtained by integration 

of compression experiment data. The smooth curve is the equilibrium entropy. The dashed line is 

the KA critical value of entropy. The compression rates are (from top to bottom) 3 ■ 10^, 10^, 10^ 

sweeps per unit of inverse System of size 20^, averages over five samples. 
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FIG. 2. Relaxation of the excess specific volume Av = v{t)—Voo after a quench to the subcritical 
value = 1/2.2. The straight line corresponds to Av oc 
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scaled time 



FIG. 3. Squared displacement R^{t + tyj,tyj) vs. scaled time t/t^ff for t^j = 10, 10^, 10^, 10^ 
and to = 10^ {t^^ = tw + to)- System of size 20^, average over five samples. 
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